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We investigate the conductivity of graphene sheet deformed over a gate. The effect of the defor- 
mation on the conductivity is twofold: The lattice distortion can be represented as pseudovector 
potential in the Dirac equation formalism, whereas the gate causes inhomogeneous density redistri- 
bution. We use the elasticity theory to find the profile of the graphene sheet and then evaluate the 
conductivity by means of the transfer matrix approach. We find that the two effects provide func- 
tionally different contributions to the conductivity. For small deformations and not too high residual 
stress the correction due to the charge redistribution dominates and leads to the enhancement of the 
conductivity. For stronger deformations, the effect of the lattice distortion becomes more important 
and eventually leads to the suppression of the conductivity. We consider homogeneous as well as 
local deformation. We also suggest that the effect of the charge redistribution can be best measured 
in a setup containing two gates, one fixing the overall charge density and another one deforming 
graphene locally. 

PACS numbers: 73.22.Pr., 73.23.Ad., 73.50.Dn, 46.70.De 



I. INTRODUCTION 

Graphene is a novel material with highly unusual elec- 
tron properties, related to the Dirac form of its energy 
spectrum at low energies, and demonstrated in many 
seminal experiements (for review seei*^) . Experiments on 
single-layer graphene have been performed on the flakes 
obtained by exfoliation as well as grown on a substrate. 

Graphene also has excellent mechanical properties. In- 
deed, the elastic properties have been measured on sus- 
pended graphene flakes mechanically deposited over a 
hole by indentation in an atomic force microscope^"—; 
the results showed that graphene is incredibly stifl[, with 
the breaking strength of the order of 40N/m, the Young 
modulus of 1 TPa, and possibility to be stretched elas- 
tically up to 20%. Bending properties have been deter- 
mined experimentally for several-layer graphene flakes^ 
and are not yet available for a monolayer. Theoretically, 
these properties have been predicted from the calcula- 
tions using the analytical form of the interatomic poten- 
tial, and from molecular dynamic studies^'^. Graphene is 
currently one of the most prospective candidates for high- 
frequency nanomechanical resonators^iS, with the quality 
factor and eigenfrequency extracted from measurements 
to be Q ~ 75, /o = 70.5 MHz for a monolayer, and 
Q « 120, /o = 42 MHz for 15nm thick graphite. Qual- 
ity factor further increases with decreasing temperature. 
An alternative method to investigate elastic properties 
of graphene is to put the film on a flexible substrate and 
deform the substrate^Siiii. The strain influences optical 
phonon spectum^^, which has been measured by Raman 
spectroscopy. 

Recent experiments combine mechanical and electrical 
properties of graphene by measuring conductivity^^"— 
of suspended graphene flakes. This is a very promising 
direction since suspended graphene flakes exhibit much 
higher mobility than graphene on substrate due to much 
weaker disorder — iH. Potentially electrons can produce 



back-action on the resonator. Graphene resonators are 
expected to have high sensitivity to mass and prebuilt 
strain^, so that they can be used to ultra-sensitive mass 
detection. 

Suspension of graphene flakes always leads to their de- 
formation, which in turn affects the conduction proper- 
ties of graphene. Deformation creates inhomogeneous 
elongation of the lattice constan t^^d'' which locally af- 
fects the electron spectrum of graphene. One way to look 
at the variations of the band structure of the strained 
graphene is to perform density functional calculations^^. 
Alternatively, the variation of the lattice constant can be 
represented at the level of Dirac equation in the form of 
pseudomagnetic fieldsi^. Ref. [13 pointed out that local 
shifts of the Fermi surface in suspended graphene in the 
vicinity of the Dirac point can block the conductivity — 
if the Fermi-surfaces at different parts of the flake do not 
overlap, the conduction is tunnel rather than metallic. 
Effects of disorder due to charged impurities and midgap 
states, optical and acoustic phonons were taking into ac- 
count for calculating conductivity of gated graphene ii*2£. 
For strong enough deformation, graphene quasiparticles 
can become localized^i. 

In experiments, graphene flakes are typically sus- 
pended over a back-gate. This gate redistributes the elec- 
tron density in the flake due to the spatial variation of the 
capacitance. The regions in the center of the suspended 
part of the flake have higher electron density then the 
regions near the clamping edges, as the central part is 
closer to the gate. This density redistribution affects the 
transmission coefficients through the entire flake. The 
corresponding effect on the piezoresistivity in ballistic 
regime is of the first order in the maximum deformation 
of the fiake in the transverse direction, and it increases 
the conductivity. This has to be contrasted with the ef- 
fect of the pseudomagnetic fields which suppress the con- 
ductivity. The contribution from pseudovector potential 
depends on the strairkii over the flake and is of the sec- 
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ond order in the maximum deformation. Thus, this con- 
tribution is expected to be weaker than effect from the 
charge redistribution. We will show however that this 
effect can be important for graphene under high enough 
residual stress. Inhomogeneous deformation of graphene 
yields the corrections to the conductivity which are of 
the fourth order in the maximum deformation, which is 
even smaller. 

In this Article, we calculate the effect of the gate- 
induced density redistribution on the conductivity of the 
graphene flake. We find that, indeed, for high residual 
stress the correction resulting from the pseudovector po- 
tential is important, and the correction to the conduc- 
tivity is negative. We mostly focus on the regime of 
low residual stress and show that the correction from the 
charge redistribution becomes the most important. 

Experimentally, influence of deformation on the con- 
ductivity would be difficult to observe on a suspended 
graphene flake with one gate since the main effect of the 
gate is the global shift of the density rather than its re- 
distribution. To separate density redistribution and elas- 
tic deformation, one needs to employ two gates. For in- 
stance, one can use the configuration with a large bottom 
gate and a narrow top gate. The bottom gate deforms the 
graphene flake and determines the maximum transverse 
deformation ^max ■ When voltage is switched on the nar- 
row top gate it does not influence much of deformation 
of the flake depleting the charge density below the top 
gate. Since the region under the top gate has the lowest 
density it determines the conductivity of the whole flake. 
If this region is brought to the Dirac point, the correction 
to conductivity is determined only by the deformation of 
graphene^'' and is proportional to {Cmax{Vg) / L)'^ , with Vg 
and L being the voltage applied to the bottom gate and 
the length of the strip under the top gate. However, for 
higher voltages the charge redistribution is more impor- 
tant, and the correction to conductivity is proportional 
to £,max{Vg)/d, d being the distance to the bottom gate. 
Instead of the top gate, one can use an AFM tip. 

The paper is organized as follows. In Section|lT]we de- 
rive equations for the deformation of suspended graphene 
from general theory of elasticity. We consider two situa- 
tions — graphene deformed homogeneously by a gate and 
graphene deformed locally by an AFM tip. The capac- 
itance between the gate and suspended graphene varies 
due to deformation of the flake. We calculate the den- 
sity redistribution over the flake taking into account the 
shape of the flake. In Section IIIIl we use these results 
to evaluate correction to the conductivity. We use the 
perturbation theory to calculate the transmission eigen- 
values, and the correction to the conductivity is obtained 
using the Landauer formula. This correction can be big 
for sufficiently strong deformations of the fiake which can 
be produced by an AFM tip. In Section IIVI we discuss 
the results and the regimes not considered in this Article. 
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FIG. 1: Side view of a deformed graphene flake suspended 
over a gate. The deformation is caused only by the interaction 
with a gate (top) or by the gate and an AFM-tip (bottom). 



II. DEFORMATION OF THE GRAPHENE 
SHEET 



In this Section, we calculate the profile of the graphene 
sheet formed by electrostatic forces induced by the gates. 
For this purpose, we decompose the total energy of the 
flake as the sum of electrostatic and elastic energies. We 
consider a graphene flake of the length L (direction x) 
and the width W (direction y) . For simplicity, we assume 
W ^ L. An undeformed sheet occupies a part of the 
plain z = 0; the electrostatically induced deflection is 
£,{x,y). Below we only consider small deformations so 
that we can stay within the limits of linear theory of 
elasticity (Hooke's law). At stronger deformations, as 
expected from the general theory^^ and also confirmed by 
theoretical modeling^i and by experiments^ on graphene, 
non-linear terms become important. However, there is a 
considerable parameter range, with the displacements up 
to 50 nm, where the linear regime is still valid. We discuss 
the terms which go beyond Hooke's law2i in Section HVl 

For electrostatic energy, similarly to Ref. [H, we model 
the system as a capacitor between the flake and the gate, 
with the distributed capacitance Cg dependent on the 
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profile of the flake, 



Cg = / v)]dxdy. 



(2.1) 



Electrostatic coupling to the leads is modelled via con- 
tact capacitances C^, Cr and resistances Rr, see. 
Fig dl]). The total electrostatic energy of the system car- 
rying the charge Q is 



CLCgVLVg 



Felectr - ^ T^TT + 7=r(CL^L + C's^s) 



2Co 



2Cn 



with Co =Cl +CR + Cg. 

From now on, we assume that the contacts are ideal, 
Cl — Cr = 0, and thus the electrostatic energy is 



F, 



electr 



(2.2) 



The effect of non-ideal contacts is discussed in Section 
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A. Elastic energy 

We evaluate the elastic energy in the thin-plate ap- 
proximation. The elastic energy consists of the bending 
contribution Fi(^(a;,j/)) and the stretching contribution 
F2{uap{x, y)), where Uap{x, y) is the deformation tensor, 
and a and /3 denote the coordinates in the plane of the 
sheet (x and y). In the linear regime, the bending contri- 
bution is less important than the stretching one, however, 
we consider both contributions for completeness. Explic- 
itly, we have^ 



FiiO 



D 
2" 



dxdy 



j j dxdy{l - 



dx^ dy^ 
\dxdy J dx'^ dy'^ 



and 



F2{Uaf}) = ho 



(2.3) 



(2.4) 



Here D = £'/ig/(12(l - ct^)) is the bending ridgity, E 
is the Young modulus, a is the Poisson ratio, is the 
thickness of the plate (graphene flake) , and ctq.^ is the 
stress tensor. 

In addition, if a local force (for instance, an AFM 
tip) acts on the graphene flake, it is best represented 
by external pressure Pext{x,y). The work of this ex- 
ternal pressure to deform the flake by d£,{x,y) is F^ = 
J PextS£,{x,y)df , where df is the surface element. 



The profile of the sheet is determined by minimizing 
its total energy. Performing the variation, we find the 
equation describing the shape of the flake. 



Peiix,y) + Pextix,y),i^.5) 



dXj3 



0,(2.6) 



with Pei(x,y) being the electrostatic pressure on the 
plate, induced by the variation of electrostatic energy 
(|2.2p . For ideal contacts, Pei{x,y) — n^{x,y)/2eo. Here 
n{x, y) is the electron density. Eq. (|2.5p is the most gen- 
eral equation for ^{x,y) in the linear approximation of 
theory elasticity. For an infinitely wide graphene flake, 
W ^ L, the deformation in the y direction is homoge- 
neous. 

At sufficiently small deformations, the tension along 
the sheet is constant over the sheet (12.61) . hoa''^^ — TSais- 
The tension T is the sum of two contributions: 



Eho 

1 -cr^ 



■AL/L 



(2.7) 



The first one, Tq, is the residual stress which results from 
the fabrication process or is induced by the ripple forma- 
tion^"^. The second contribution, Th, is an internal force 
due to the relative elongation AL/L (Hooke's law). If 
we take this term into account, we can go beyond the 
thin-plate approximation and consider deformations big- 
ger then the thickness of the graphene layer. 

In the two following Subsections, we solve the above 
equations for two specific situations: homogeneous ex- 
ternal force (which can be produced by a bulk bottom 
gate), and local force (produced for example by an AFM 
tip). 



B. Homogeneous force: Deformation by a bottom 

gate 

Applying a voltage on a bottom gate is a standard way 
to vary electron density in graphene. If the suspended 
graphene fiake is charged, it is subject to a mechanical 
force proportional to the charge density. If the area of 
the gate is much larger than the area of the fiake, the 
electron density induced by the gate is constant almost 
everywhere, n — Q/WL, except for the clamping points 
of the flake, where it is determined not only by the so- 
lution of the Poisson equation (providing singularities at 
the capacitor edges), but also by the metallic leads to 
which the flake is clamped. Indeed, experimental evi- 
dence for this charge inhomogeneity exist and can be ac- 
cessed by asymmetry of the Dirac peak in conductivity 
However, these density inhomogeneities at the clamping 
areas very little affect the deformation, since the displace- 
ment vanishes at the edges of the flake. Therefore we can 
approximate the effect of the gate by homogeneous elec- 
trostatic pressure over the flake, P = eoV^/2(i^, Vg and d 
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being the gate voltage and the distance to the gate. The 
profile of the graphene sheet is found from the equation 



(2.8) 



where the stress T is constant over the sheet (|2.7p and 
the deformation-dependent contribution to it depending 
has to be found self-consistently, 



Th = 



Eho 



2(1 



(2.9) 



(the case for inhomogeneous Th derived in Ref. [2J is 
discussed in Section IIVI and does not induce significant 
difference in results). For the boundary condition cor- 
responding to the clamping the sheet, ^(0) — ^{L) — 
0,0) = C{L) = 0, the profile is 



sinh fiL 




(cosh iix — I) ~ sinh fix + jix 



(2.10) 



The profile (|2.10l) is parabolic in the middle of the strip 



(as noted in Ref. |17D . As we show below, in graphene the 
dimensionless parameter fj,L assumes large values. In this 
case, the profile can be simplified, and near the middle 
of the strip has the form 



e(x) 



PL 
2T 



X 

T 



(2.11) 



Close to the edges, the profile becomes ^{x) — 
PfiLx^/iT. Substituting this shape into Eqs. (|2.3p and 
(|2.4p . we find the values of the parameters Fi and F2, 



and 



Fi = 



P^LW 



16T^2(-8 + /iL)^ 



T^WL{l-a^) 

Ehn 



(2.12) 



(2.13) 



The maximum vertical displacement obeys the equation 

PL2 



8iTo+8Ehoemajm-cT^)L')y 



(2.14) 



The deformation of the sheet leads to the redistribu- 
tion of the electron density, which in the Thomas-Fermi 
approximation is n{x) = Vgeo/{d — ^{x)). In its turn, 
the density redistribution affects the profile of the sheet, 
and needs, in principle, to be calculated self-consistently. 
However, as soon as the displacement Cmax is much 
smaller than the distance to the gate, the later effect is 
insignificant (of the order of ^ax/d), and we will use the 
shape (|2.10p not modified by the density redistribution. 




FIG. 2: Dependence of the mEiximum deviation on the 
gate voltage, S^maxiVg). The solid curves represent the self- 
consistent solution of nonlinear coupled equations for the de- 
formation of the flake and the charge induced by the gate, 
Eqs. (|2TT4|) and (pTTe]) . The distance to the gate is d = 300 nm 
(top panel) and d = 100 nm (bottom panel). Other parame- 
ters of the graphene flake, length L = 1 fim, Young's modulus 
E = 1 TPa, Poisson ratio cr = 0.15, and the thickness of the 
flake hn = 0.34 nm, are chosen in order to model real exper- 
imental data. The resuhs are given for the different values 
of the residual stress: the curve 1 is for To — 0.001 N/m, 
2 is for To = 0.01 N/m, 3 is for Tq = 0.1 N/m. For each 
value of the residual stress, the asymptotic curves at low gate 
voltages (I2.17[) are shown as dashed lines, the curve 3 for 
high residual stress coincide perfectly with its assymptote. 
The asymptotic curves for low residual stress, Eq. (|2.19|) . are 
shown by dashed-dotted lines. The correspondence between 
the solution of equations and asymptotics for low residual 
stress is not perfect. The reason is that the asymptotics 
are calculated for the linear charge-voltage dependence, and 
n(Vg) is non-linear according to Eq. (|2.16l) for sufficiently high 
gate voltages on the flake. 



The charge over the graphene flake is determined from 
minimization of the total energy of the system with re- 
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3d \ 2 J dfiL 3 d 

(2.15) 

where the maximum deformation of the sheet in the mid- 
dle, ^max (|2.14p , and depends on charge density. The sec- 
ond term in the brackets and 1 from the third term come 
from electrostatic energy and originate from the redistri- 
bution of the charge density due to variation of the dis- 
tance between parts of deformed graphene and the gate. 
The rest (1/2) of the third term comes from bending en- 
ergy. The fourth term takes into account dependence of 
the stretching force T over the flake on the charge den- 
sity via the deflection (Eq. (j2.9p ). Calculations are made 
under the assumption /iL 3> 1, which is realistic for avail- 
able experiments. Simplifying Eq. (|2.15p . we obtain 



(2.16) 



At low gate voltages, Eq. p.l6p yields the linear gate 
voltage dependence of the electron density, uq = Vgeo/d. 
There is non-linear deviation from this dependence at 
higher gate voltages and at rather low initial strain Tq. 

The maximum deformation can be expressed analyti- 
cally in two limiting cases. First, if the residual stress Tq 
is stronger than the induced stress Th , it mostly accounts 
for the deformation of the sheet, 



^2 ^ ^0: ' 



24(1 - a2)T2 



V24(l -cr2 



Th 

^max 
n — Hq 

no 

In the case of low residual stress, one obtains 



1 fEhoP^L^'^ 
To « - - ' ° 




«To, 
(2.17) 
(2.18) 
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, (2.19) 
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_12j5^4i^. (2.20) 
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The maximum deviation ^max, obtained from the 
numerical solution of coupled nonlinear equations 
Eqs. (|2.14p . (|2.16p . as well as asymptotic expressions 
(|2.17|) and (|2.19l) . are shown in Fig. [5] for different values 
of initial stress Tq^. According to Eqs. ([2T8)) and (g^Ol, 
the nonlinear part of the charge induced on the graphene 
flake follows the dependence ^maxiV)/d- Consequently, 
we encounter several regimes for the deformation. 



large 

A 



regime 



c 
o 

c 

S Linear 
^ charge 



/ 



small -^0 



gate voltage, F 



An exact solution 
of electrostatical 
problem 

f /large 

deformation 

nonlinear cliarge 
reg ime 



FIG. 3: Regimes of the deformation of suspended graphene. 
For large residual stress To the asymptotics Eq. 1)2. 17|) are 
valid, and the charge on the flake follows the gate voltage lin- 
early. For small To the asymptotics Eq. 1)2. 19|) are valid, the 
charge is linear with the gate voltage at low gate voltages and 
starts to follow non-linear dependence with increasing voltage. 
At high gate voltages, when the deformation of the flake is 
bigger than S^max/d ~ 0.1, one needs to solve self-consistently 
the electrostatic problem and the problem of elasticity, ana- 
lytical results for this region are not available. 



• at large Tq the charge-voltage dependence is lin- 
ear for realistic parameters because the maximum 
deformation is not too large for realistic character- 
istics of graphene flake. It is shown in Fig. [5] for 
To = 0.1 N/m that the maximum deformation is in 
a good agreement with Eq. (|2.17l) : 

• at small To and low gate voltages Vg the charge- 
voltage dependence can be in linear regime, and the 
maximum deformation follows Eq. (|2.19p . We il- 
lustrate this for the flake with the parameters To — 
0.001 N/m and distance to the gate d = 300 nm 
(See. Fig. [21 top), where the solutions of cou- 
pled electrostatic and elastic equations, Eq. (|2.14p . 
(|2.16p . follow asymptotic expression Eq. (|2.19p . 
The charge redistribution does not need to be taken 
into account; 

• at small To and large Vg the system is in the non- 
linear charge regime. This situation can be real- 
ized for small distances to the gate when the cou- 
pling of the graphene sheet to the gate is large, so 
that it is possible to create large deformations us- 
ing low gate voltages. For example, at d = 100 nm 
the non-linear charge regime influences the defor- 
mation already at voltage = 2 V, at the bot- 
tom plot Fig. [2] we can see the intersection of the 
asymptotical curve (|2.19p and the actual solution 
of Eqs. dsn, (PIB . 

The schematic representation of these regimes is shown 
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in Fig. H 



C. Local force: Deformation by an AFM tip 

Next, we consider a concentrated force acting on 
graphene. This force can be provided, for example, by 
an AFM tip. The effect of the tip is modeled by strong 
pressure exerted on a narrow area of the width Z <C i. 
We assume that the problem is still homogeneous in the 
y-direction, what simplifies the calculations enormously. 
Inclusion of a pressure action in a narrow circle, which is 
experimentally relevant for an AFM tip, is not expected 
to bring qualitatively new features. We consider pressure 
P(x) = P2,0 < X < L/2 - 1/2, L/2 + 1/2 < X < L and 
p\x) = P\,L/2 -1/2 < X < L/2 + 1/2. Here P2 is the 
local pressure, and Pi <^ P2 can describe homogenious 
pressure due to electrostatics. 

The maximum displacement of the flake (realized at 
the central point) is easy to write down for /iL ^ 1 and 



Pi((Aii/2)2e^'/2-2e-^') 



P2{e 



li'^lL/A) 



IJL^T 

For Fi < P2 and 1 < e'^'-^'^Lnlfi/A we obtain 

P2IL 



4T 



(2.21) 



(2.22) 



The profile of the graphene sheet in this approximation 
becomes 



L 



(2.23) 



In the limits of weak and strong residual stress the 
deformation is determined by 



8(1 - ^2) V To 
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P2IL 



(2.24) 
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1 (EhoPiP 

2 V 1 - 
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1 fEhoPiP 
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^2^(1 



8Eho 



1/3 



(2.25) 



The dependence of the maximum deformation on the 
applied external local force is shown in Fig. The de- 
formation produced by this force is much bigger than the 
deformation caused by electrostatic pressure of the gate. 
The electrostatic problem for this case can be solved sep- 
arately from the problem of elasticity. 




FIG. 4: Dependence of the maximum deformation ^max on 
the applied force P2I for the case of the point force in the 
middle of the graphene sheet. Only the curve for To = is 
shown, since the residual stress is not important for this case: 
The strain created by deformation becomes large (more than 
0.1 N/m) already at moderate deformations in the middle, 
^max ~ 10 nm. Other parameters of the flake axe—. 



III. PIEZOCONDUCTIVITY OF GRAPHENE 
FLAKE 

It was shown experimentallyl^ that suspended 
graphene flakes are described with good precision as 
purely ballistic. Theoretically, conductance is deter- 
mined by Landauer formula^^ 



4eV/i J2 



(3.1) 



where Tn is the transmission eigenvalue in the trans- 
port channel n, and the factor 4e^//i. is conductivity of 
a single transport channel which takes into account val- 
ley and spin degeracy. The number of open transport 
channels N — Wkp/ir is proportional to the Fermi mo- 
mentum kp — {TTn/ e)^/"^ , and thus the conductivity is 
proportional to the square root of the electron density n, 
a (X \fn. 

The conductivity of graphene flake suspended over a 
gate can deviate from this dependence. To start with, 
due to electrostatic interaction with the gate, the den- 
sity becomes inhomogeneous^^. In particular, Poisson 
equation leads to the square root divergence of the elec- 
tron density at the clamping points, as in any capacitor 
(see e.g. Ref. HI). To treat this divergence properly, one 
has to take into account electrostatic interaction with 
the contacts near the edge of the graphene strip, which 
modifies significantly the electron density near the edge, 
removing the divergence. However, the effect of this in- 
homogeneous density close to the contacts does not affect 
the piezoconductivity of the flake, since the deformation 
close to the clamping points is very weak, and thus it can 
be included into the contact resistance at the clamping 
points. 
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We now turn to the effects of the deformation on the 
conductivity. Deformation of graphene can change the 
conductivity by inducing changes in the band structure 
(which resuhs in pseudo-magnetic fields) as weh as by 
changing the electron density over the fiake. We consider 
both these mechanisms and will show that typically the 
effect of the density redistribution dominates. 

Electrons in graphene obey Dirac equation. Deforma- 
tion of the fiake infiuences on the Dirac equation in three 
ways — it shifts the K-points by a certain amount Sk/kp 
(pseudomagnetic field), renormalizes the Fermi velocity 
by Svp/vp, and induces the variation of the electron den- 
sity on the fiake 5n/n. The deformation correction to the 
conductivity is thus a function of these three dimension- 
less parameters. 

The pseudomagnetic field, produced by the shift of the 
/C-point, is caused by stretching and bending. The shift 
of the K-po'mt due to stretching generates the vector 
potential^^'^^ 
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(3.2) 
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FIG. 5: Schematic behavior of piezoconductivity. For low 
residual stress To the correction is mainly due to the charge 
redistribution and has positive sign. For high residual stress 
the correction is negative. 



where C is the order of 1, and /3 — — 91og(i)/91og(a), 
t and a being the overlap integral in the tight-binding 
model and the lattice parameter, respectively. For L <^ 
W one has Uxy = 0, and hence A^*^ — 0. The defor- 
mation is homogeneous within the limits of applicability 
of Hooke's law, and thus u^x = const and Ay^^ = const. 
This means that there is no pseudomagnetic field over the 
graphene flake. The pseudomagnetic field only appears 
in the region where the fiake goes from the substrate 
to the suspended stateii and, as noted above, its effect 
to the piezoconductivity is small, of the second order in 
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(3.3) 



where the deformation on the edges has been estimated as 
Uxx = ^^ax/L^+TQil~cr^)/Eho, and kp = ^7reoVg/de. 
Taking into account the value of C[3/ci^, we obtain 



SCTK 



= 205, 



r2 



+ 



To(l 



Eho 



(3.4) 



Note the contribution from two terms induced by defor- 
mation stress and residual stress, as well as multiplication 
with the big prefactor 205. 

The underlying physical picture for the model of 
Ref. |T3 is that the graphene fiake is " glued" to the walls 
at the suspension point. Whereas this has been realized 
in some experiments it describes the situation when the 
residual strain Tq is of the same order or higher than the 
strain induced by the gate voltage. The residual strain 
results from the fabrication process and is most likely to 
be created by impurities in the substrate. It can be made 
low on purpose since the strain is reduced after anneal- 
ing^. In the opposite situation, when the residual stress 



is not significant, the pseudomagnetic field is inhomoge- 
neous and distributed over the whole suspension area. 

The pseudomagnetic field is also inhomogeneous if one 
considers the bending contribution. Bending leads to 
the inhomogeneous modification of the overlap of the or- 
bitals, and the resulting pseudovector potential has the 



bend 



A ben 

Ay 



tbend fO^{a,x) e%^,x) 



(3.5) 



6{a,x) being the angle between normal vectors to the 
graphene surface at the points x and x + a, and the con- 
stant-^^ thend = 3.21—. The shape dependence of 9{a, x) 
has the form 



^(a,a;) 




This yields A. 



bend 



my 



{^tbend/S>a){S,maxa/L^r. Hcncc the 



contribution from bending is approximately [a/ L) times 
smaller then from stretching without residual stress, and 
is thus negligibly small, even though the resulting mag- 
netic field is not homogeneous. 

The easiest way to estimate inhomogeneous stretch- 
ing of graphene is to take Hooke's law in the local form, 
Th{x) — EhoUxxix). Since the maximum relative de- 
formation can be estimated as Uxx = Cmaa:/^^' naively, 
the correction from non-homogeneous stretching is of the 
same order as the one from delta-functional pseudomag- 
netic field at the clamping edges. We show below, how- 
ever, that the correction from non-uniform stretching is 
of the order of (.tiax/^^^ but still due to large prefactor 
it can reduce the conductivity at low gate voltages. 

Another effect induced by the deformation is the renor- 
malization of the Fermi velocity. The renormalized 
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FIG. 6: The correction to conductivity. Parameters of the 
flake are the same as for Fig. [2] (top). Asymptotic expres- 
sions for the high gate voltage are shown by the dashed-dotted 
lines. Additionally, the correction due to the delta-functional 
pseudomagnetic field at the suspension regionsi^, Eq. (|3.4[) . 
is shown (with the opposite sign). 



value of the velocity can be derived from the tight- 
binding model. Assuming that the graphene sheet is 
only deformed in the x-direction, we find that the x- 
component of the Fermi velocity is unchanged wheres 
the y-componcnt is renormalized, 



Vpy =Vp{l- Cjiu^.^), 



(3.6) 



so that approximately vpy ~ wf(1 — ■fmax/-^^)- The effect 
of the renormalization on the conductivity is not signif- 
icant and has the order of magnitude ^max/^'^- Note 
that this is the same dependence on £,max/ L as for pseu- 
domagnetic fields, however, it is not enhanced by a big 
prefactor. 

The influence on conductivity of such change in fermi- 
velocity is not significant. This influence can be in 
principe measured experimentally as the conductivity 
variation at the Dirac point, similarly to how we explain 
below in Subsection IIIIBI 



Correction to conductivity due to the charge 
redistribution 



Redistribution of electric charge due to interactions 
with the gate is found from the assumption that the 
potential along the graphene sheet is constant, U{x) = 
5Q{x) 1 5C{x) = const, where 5C{x) is the capacitance 
of the element of the length 5x of graphene, 5C{x) — 
WSx/4:n{d — £,{x)), and SQ{x) = n{x)WSx is the charge 
of this element. In the first order approximation, this 
gives dn{x)/no = £,{x)/d. 

The conductivity of graphene is proportional to charge 
density n, and thus the contribution to conductivity 
due to charge redistribution is expected to be linear in 
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FIG. 7: Dependence of the piezocorrection to conductivity on 
the gate voltage for fixed deformation, obtained by solving the 
Dirac equation by exact transfer matrix method. Here, Sau, 
the curve marked 1, is the correction with only charge redis- 
tribution taken into account, and SaA,u, the curve 2, encom- 
passes both contributions, the one due to non-uniform tension 
and the one due to charge redistribution. At low gate voltages 
the correction S(ja,u is mostly caused by pseudo-magnetic 
field and is negative, for higher gate voltages it changes sign 
and approaches to 5au- Parameters of the graphene flake 



the maximum deviation from the homogeneous density, 
Srimax- Thus, the correction to conductivity is expected 
to be Sa/a ~ ^max/d. Before starting the calculation 
of the correction to conductivity, we estimate the range 
where this correction of the order of S,max/d is more im- 
portant than the correction due to pseudovector potential 
which we considered above. Fig. [Sj 



ro[N/m] < 10-3^ 




(3.7) 



As noticed in Ref. [iTl the pseudovector potential at 
low gate voltages blocks conductivity, this is seen from 
Eq. (|3.4p . For large deformation the expressions (|2.19p 
are valid, residual stress is not important any more, and 
thus the gate voltage should be large enough to see the 
decrease of conductivity. 



K,[V] >2.i 



1 



d[/im] 



(3.8) 



This deformation is so strong that in can not be reached 
in practice. 

For the deformation with AFM the residual stress is 
not important and at deviations 



Cmaa;[nm] > 2.5 



d[/im] y d[/im] 



(3.9) 



correction to pseudovector potential starts to suppress 
the conductivity. 
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To calculate the correction due to the charge redistri- 
bution, we notice that the density variation is translated 
into the correction for conductivity via the variation of 
the transmission probabilities T„ , which are the eigenval- 
ues of the matrix iH, i being the transmission matrix of 
the graphene sheet. The transmission eigenvalues tq are 
determined in Appendix by the transfer matrix method. 
The correction to the conductivity is linear in the density 
shift Sn, and consequently in the maximum deformation 
^max (as is shown above from simple qualitave consider- 
ations). It has the form (see Appendix) 



Sau = ^4|i,|^^sinfcL 



where 



X / dx^^^smk(L- x)smkx, (3.10) 
Jo d 

(cos^ kL+kp sin^ kL/k^)~^ is the transmis- 



sion probability for the mode labeled by the transverse 
momentum q = 2T:n/W, n being an integer number, and 
fc is a wave number in the direction along the strip, so 
that kj, = k'^ +q^. 

To carry out more detailed analysis, we consider spe- 
cific deformation setups discussed in Section — homo- 
geneous and local deformation. 

Eq. p.lOp can be analyzed analytically for small and 
large values of the parameter kpL, which characterizes 
the charge density over the flake. The correction to 
conductivity for the homogeneous deformation (bottom 
gate) has the following asymptotic behavior for small and 
large values of kpL (for more details, see Appendix), 




1 < kpL; 
.{kpLf/d, kFL<^l. 



(3.11) 



Taking into account the functional dependence of the 
maximum deviation for small and large initial stress Tq, 
Eqs. (|2.19p and (I2.17p . we get the asymptotic dependence 
of the correction to conductivity on the gate voltage, for 
Th < To: 




L^V^/Tod\ 



1 < kpL] 
kpL < 1, 



^4/3^2/3/^5/3^ l^kpL- 

{VL^/d^f'\ kpL^l. 



(3.12) 



(3.13) 



Fig.|6]shows the exact result of summation over modes 
Eq. (|3T0| . At both high {kpL > 1) and low {kpL < 
1) gate voltages, the correction follows the asymptotic 
behavior both for weak and strong residual stress Tq, 
Eqs. (|3.13l) and p.l2|) . On the same plot we compare 
the correction we found with the correction due to pseu- 
domagnetic fields at the edges^. The latter one has a 
different sign (conductivity decreases with the stress). 
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FIG. 8: Dependence of the piezocorrection to conductivity on 
the maximum deviation, 5a{(^max) / o , at the fixed gate volt- 
age Vg = 0.03 V (top), 3 V (bottom). We show both the 
correction related to the term due to non-uniform pseudovec- 
tor potential, and the contribution without this term. For 
Vg = 0.03 V, we also include the best HtV — a^max + l^^max 
which represent the sum of Hnear in £^max correction due to 
charge redistribution and the correction due to nonuniform 
pseudovector potential P^f^^^. Parameters of the graphene 
fiake are^-. 



For high residual stress this correction is more important 
than due to charge redistribution, according to Ref. [l3 
it can block conductivity. For low residual stress it is 
about 10 times lower than the increasing conductivity 
correction. The oscillations of da /a have the period of 
kpL d are associated with the shift of Fabry-Perot res- 
onances in conductivity for deformed graphene flake as 
conpared with an undeformed flake. This shift occurs 
since the effective longitudinal wave vector of an electron 
in graphene depends on the deformation since it feels 
different charge density over the graphene flake. Note 
also that the contribution from pseudomagnetic fields 
does not oscillate since the value of kp is the same for 
the whole flake. The flrst order perturbation theory in 
S,max/d is valid until this parameter reaches a rather large 
value, ^max/d ~ 0.1 (see Appendix for more details). 

For the case of local deformation, using the graphene 
profile (|2.23p and using the same technique as in Ap- 
pendix, we find the correction for conductivity due to 
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the charge redistribution, 



6a 



1 < kpL; 



a I QmS^r,^ax{kFLf/d, kpL < 1. 



(3.14) 



Note that the asymptotic behavior for large kpL has 
the sarae form as for homogeneous deformation. 

We can also estimate the influence of inhomogeneous 
pseudomagnetic field assuming the local form of Hooke's 
law as in^^ and using the perturbation theory for the 
transfer matrix, as detailed in Appendix. We find that 
the first order perturbation theory correction in pseu- 
dovector potential vanishes, whereas the second order 
correction can be estimated as 



5a I 



5.5 • 10' 



L4 



(3.15) 



At low gate voltages and large deformations (for instance, 
induced by local deformation), this correction can be 
more important that the one from the charge redistribu- 
tion, and thus the conductivity will be suppressed. From 
comparison of Eqs. (|3.15l) and (|3.14l) this supression hap- 
pens for deformations: 



L 



1/2 



(3.16) 



We demonstrate this by solving numerically by transfer 
matrix method the Dirac equation with additional poten- 
tial due to charge redistribution and pseudovector poten- 
tial. Fig. [71 At fixed large £,max — 12 nm (estimated using 
Eq. p.lGp ) and at low voltages the conductivity starts to 
decrease due to inhomogeneous tension distribution in 
the flake, and at higher voltages increases again due to 
the effect of charge redistribution. Fig. |S] shows that for 
small gate voltages lower maximum deformation S.max is 
required to reach the point where the conductivity starts 
to decrease, in agreement with Eq. p.l5p . For high gate 
voltages Vg ^ 3 V pseudomagnetic fields lead to satura- 
tion of the conductivity rather than to its decrease. 



B. Two-gate geometry 

Conductivity can also be used to measure relative 
stretching of deformed suspended graphene. Note 
that the influence of stretching on the conductivity of 
graphene deposited on a substrate has been demon- 
strated experimentally^'^. For suspended graphene it is 
more difficult to extract the value of stretching than 
from the graphene on the substrate, since the gate volt- 
age simultaneously varies the concentration and deforms 
graphene, as shown above. 

To measure relative stretching of suspended graphene, 
we propose the two-gate geometry (Fig.[S]). The deforma- 
tion of the graphene flake is created by the large bottom 
gate, the influence of the top gate on the deformation is 




FIG. 9: Schematic dependence of conductivity of suspended 
graphene on the top gate voltage for several fixed bottom 
gate voltages. The conductivity at the Dirac point is slightly 
shifted due to change in Fermi velocity caused by deforma- 
tion. The difference between the values of conductivity of 
the Dirac peaks for different bottom gate voltages, Vbi and 
142, is proportional to the difference in relative deformations, 

(o-o(Vbi) - CrD(V(,2))/o"D ~ Uxx(Vbl) - Uxx{Vb2)- 



small as the top gate is narrow. The top gate is used to 
vary the charge density in the region underneath it. In 
this geometry at the fixed voltage at the bottom gate one 
can move through the Dirac point by varying the volt- 
age at the top gate (the experiment for bilayer with two 
gates on the substrate'^'*). The value of conductivity at 
this point depends on the deformation. 

The stretching of the graphene flake, as discussed 
above, induces variations of the conductivity for two rea- 
sons. First, it induces pseudomagnetic fields. These, 
however, can be gauged away of Dirac equation'^^ at the 
Dirac point and do not influence the conductivity. Sec- 
ond, it shifts the Fermi velocity. The relative shift is pro- 
portional to the deformation, Svp/vp ''^ ^max/^'^^ ^^'^ 
leads to the positive correction of the conductivity at 
the Dirac point, Sa/a ^ Svp/vp ~ ^max/^'^- Thus, 
for different bottom gate voltages, which is equivalent 
to different maximum deformations S,max, the conduc- 
tivity at the Dirac point is slightly different, and the 
relative graphene stretching can be restored from this 
dependence. For example, consider the dependence of 
conductivity on the top gate voltage for different bot- 
tom gate voltages (Fig. [S]). At a flxed value of the bot- 
tom gate voltage, the conductivity as a function of the 
top gate voltage exhibits a peak dependence, with the 
minimum corresponding to the Dirac point. The differ- 
ence between the values of conductivity at Dirac peaks, 
an, for different bottom gate voltages, Vti and Vb2, is 
proportional to the difference in relative deformation, 
(CTu(Vbi) - aD{Vb2))/aD ^ Uxx{Vbi) - Uxx{Vb2) (we re- 
member that Uxx ~ ^max/-^'^)- 
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IV. DISCUSSION 

In this Article, we investigated two mechanisms which 
affect the conductivity of suspended graphene — charge 
redistribution induced by the gate(s), and pseudomag- 
netic fields induced by the deformation of graphene. We 
find that for the small residual stress Tq, the charge re- 
distribution mechanism dominates. For low gate voltages 
and strong deformation, which experimentally is best re- 
alized by using AFM, the correction due to nonuniform 
pseudomagnetic fields is more significant. The correction 
due to pseudovector potential at the region of suspension 
can decrease conductivity at the large residual stressA^. 
It is important that the two mechanisms provide correc- 
tions to conductivity which are of different signs. Indeed, 
the correction from pseudomagnetic fields suppresses the 
conductivity^^ by shifting K-points due to the vector po- 
tential. The shift is different at different points of the 
suspended sample, and if the deformation is big enough, 
the Fermi circles at the clamping points and at the centre 
of the flake do not overlap: The system becomes insulat- 
ing. If now we take into account the effects of the gate, 
not only the Fermi circles are shifted, but their radii are 
greater at the center of the flake since the charge density 
is greater in the areas closer to the gate. The increase of 
the radii and the shift of the center compete, and we find 
that typically the radius increase is more important. 

It is difficult to measure piezoconductivity only by us- 
ing a bottom gate since the gate voltage not only bends 
graphene and produces the correction to the conductiv- 
ity, but also shifts the overall charge density. The density 
dependence of the conductivity is different from the den- 
sity dependence of the correction. Thus, to extract the 
value of piezoconductivity, one has to compare the con- 
ductance of deformed and undeformed graphene sheet at 
the same density, which can only be done in the one- 
gate geometry by comparing the results with the theo- 
retical prediction. In contrast, the two-gate setup, with 
a bottom gate fixing the overall density and the top gate 
(which can be an AFM tip) creating the deformation is 
more convenient to extract piezoconductivity. One can 
fix the voltage on the bottom gate and start to deform 
the fiake with the AFM tip. At low gate voltages the 
conductivity decreases due to the pseudomagnetic fields, 
whereas at higher voltages it starts to grow due to the 
charge redistribution. 

In the real experimental situation, the AFM tip has a 
point shape, whereas in this Article we considered for il- 
lustration the deformation homogeneous in one direction, 
i.e. replaced the tip by a rod. Non- homogeneous defor- 
mation in all directions creates pseudomagnetic fields, 
with the conductivity depending not only on the trans- 
verse displacement, but also locally on the position over 
the graphene sheet. The conductivity is the largest if the 
tip is placed in the middle of the sheet, and decreases if 
the tip moves to the side. We can understand this behav- 
ior from a simple reasoning. Indeed, the electrons which 
from the two sides of the tip feel the pseudomagnetic 



fields and interfere similarly to an Aharonov-Bohm ring. 
The interference is more destructive if the tip is further 
from the center, and thus the conductivity decreases. 

Another parameter which affects the conductivity is 
the residual stress Tq. It can be varied experimentally 
for instance if one uses graphene suspended over piezo- 
substrate. Putting voltage on the substrate would induce 
extra stress on graphene, and one can move from the situ- 
ation where pseudovector potential blocks the conductiv- 
ity at low gate voltages to the case where residual stress 
does not play a role and the correction due to charge 
redistribution increases the conducitivity. 

In this Article, we considered ideal ballistic graphene. 
In particular, we disregarded the contact resistance, as- 
suming the clamping points to be ideal contacts. Finite 
transparency of the contacts would suppress both the 
conductivity itself and the piezocorrection to the con- 
ductivity; in addition, it would raise the amplitude of 
Fabry-Perot resonances. 

For strong deformations of the graphene sheet, the 
problem becomes much more complicated, since one has 
now to solve elasticity equations self-consistently, tak- 
ing into account that the displacement depends on the 
charge redistribution. This leads to additional terms in 
the equations of the elasticity theory. Taking into ac- 
count infiuence of the density redistribution on the term 
with electrostatic pressure in the equation of deforma- 
tion, one can show that the self-consistency condition 
increases the deformation in the middle of the graphene 
sheet. This effect only becomes important at sufficiently 
strong deformations. 

Finally, we assumed that undeformed graphene is flat. 
In reality, it is always rippled, and, in principle, one needs 
to use the elasticity theory for membranes. However, we 
do not expect that taking ripples into account would sig- 
nificantly affect the results of this paper. First, the rip- 
ples are small and have a large radius of curvature, which 
means they are very little affected by the overall defor- 
mation of the graphene sheet. Second, the main effect of 
the ripples is to renormalize the energy over the graphene 
sheet^. We thus expect that our results are valid, but 
for renormalized energy over the flake (energy is deter- 
mined by gate voltage in clean case, and is renormalized 
in the rippled case). 



V. ACKNOWLEDGEMENTS 

We thank A. Akhmerov, C. Beenakker, M. Fogler, S. 
Goossens, F. Guinea, A. Isacsson, and A. Kharimkhod- 
jaev for discussions. We acknowledge the financial sup- 
port of the Future and Emerging Technologies program 
of the European Commission, under the FET-Open 
project QNEMS (233992), of the Dutch Science Foun- 
dation NWO/FOM and of the Eurocores program Euro- 
Graphene. 



12 




FIG. 10: Dependence of the relative correction to conductiv- 
ity on kpL, 5cj / oikpL) for constant ^max/d = 1/3000, for the 
correction of the first order in ^max/d, Eq. (|6.6() . the curve 1, 
and exact transfer matrix solution of integral equation, the 
curve 2. The correction from the exact solution has the same 
dependence on kpL as the first order correction, the oscilla- 
tions are in the same phase. Asymptotics for small and large 
kpL, Eq. p.lOf) . are shown as dashed lines. Parameters of 
the fiake are^. 



VI. APPENDIX. PERTURBATIVE 
CORRECTIONS TO CONDUCTIVITY 

In this Appendix, we calculate the corrections to the 
conductivity due to both charge redistribution and pseu- 
domagnetic fields, using the preturbation theory. 

The Dirac equation for one valley in graphene has the 
form 

VF^P + SU{x,y) = E , (6.1) 
with a = {(7x,cry), p = ipx,Py), 

Px = -ifldx + Ax,Py = -ihdy + Ay, 

Ax{x,y) and Ay{x,y) being the components of the pseu- 
domagnetic vector-potential^, given by Eq. (|3.2p . and 
SU{x,y) is the additional electrostatic potential due to 
the charge redistribution over the graphene flake. It 
is determined by local variations of the Fermi energy 
over the flake. Since the Fermi energy depends on 
the charge density over the flake, Ep{x) = Tivpkp{x), 
kp{x) = yjiin^x) / one has 

5U{x)/E = 6kp{x)/kp = 5n{x)/2n = ^{x)/2d. 

We only consider the deformation homogeneous in y- 
direction.Then both Ay and dU only depend on the co- 
ordinate X, and Ax = (see Section UlI)) . The problem 
becomes effectively one-dimensional since the momentum 
q in y-direction is conserved. It is convenient to use the 
transfer matrix representation of Dirac equatior*^ to cal- 
culate the correction to the conductivity caused by the 



deformation Ay{x), SU(x), 
Th {x2,xi) =%h{x2,xi) - (6.2) 
dx%H{x2,x) {(jJU(x) +iaxAy{x))%i{x,xi) , 



where Th is the Hadamard transformed transfer matrix, 
and Tqh is the Hadamard transformed transfer matrix of 
the unperturbed system, 



%H = exp {ia^kpL + a^qL) . 



(6.3) 



We perform the perturbation expansion of the integral 
form for Eq. (|6.3|) . and in the first order in dU{x) and 
Ay{x) we obtain 

Tl{x2,Xi) = To{x2,xi) -(6.4) 



dxToix2,x) {az6U{x) + iaxAy{x)) %{x, xi). 



The conductance of the graphene sheet is determined 
by Landauer formula (|3.1|) . According to general scat- 
tering theoryi^'', the transmission matrix element t is an 
inverse element of Th, 



(6.5) 



Taking into account Eq. (16. 5p . Landauer formula p.ip . 
and the definition (|6.5p . the first order corrections to con- 
ductivity due to electrostatics and pseudo-magnetic field 
are 



■ sin fcL, 



lu 



— sin k{L — x) sin kx^ 
L/ Zcl 



(6.6) 



(6.7) 



I A = 



dxSA{x) X 



X (sin^ fci — 2cosfcLsinfca;sinfc(L — a;)), (6.8) 

where q — 2T:n/W is a wave vector in the y-direction, n 
is an integer number, and fc is a wave vector along the 
strip. 



k^+q^ 



kp. 



Furthermore, tq is the transmission probability for clean 
system for the mode g, and 

\tq\^ = (cos^ kL -I- k\ sin^ kLjk^Y^. 

Note that the first-order correction due to the pseudo- 
vector potential (16. 7p only contains odd powers of g, so 
that the sum over q vanishes. Thus, the first-order cor- 
rection to the conductivity is determined solely by the 
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density redistribution. It is linear in the maximum devi- 
ation S^rnax/d for Small deviations. 

First, we remark on the validity of Eq. (16. 6p . The 
expansion of the expression 



1 - 'i{tqt\YlukFL^^-^ sin kL 

has been made under assumption that the second term is 
small in comparison with unity due to the small prefactor 
£.maxkFL/d. Following this argument, the expression for 
the first order correction to the conductivity in ^max/d is 
formally only valid for S,maxl dkpL <§C 1. However, solv- 
ing the integral equation numerically, we find that this 
expression is valid for a broader parameter range. We 
compare results of calculations for the first order cor- 
rection Eq. (|6.6p and numerical solution of Eq. (|6.3p for 
the two cases: for the fixed ratio ^max/d and for the 
fixed value of kpL. For the first case, the dependence 
of Serf a on kpL shows the same oscillation period and 
the same asymptotic behavior at large kpL, Fig.[TOl For 
the second case, at large kpL ~ 40 (for the distance 
to the gate d = 300 nm this corresponds to the gate 
voltage = 3 V) the expansion clearly ceases to be 
valid, see Fig. [TTJ We thus conclude from the results of 
our numerical solution that the expression for the cor- 
rection linear in (,max/d is applicable until ^max/d <C 1, 
which is weaker than the perturbation theory suggestion 

^maxkpL/d-^ 1. 

The second order correction to conductivity contains 
also a term with the pseudo- vector potential, the magni- 
tude of the term being {£,max/L)'^. 

We consider both corrections separately. Now we per- 
form the analysis of Eq. (|6.6|) for deformation with con- 
stant pressure. For this case, the shape of the strip is 
nearly parabolic (Section |Tl| and can be approximated 
as 



'-{x-L/2)^. 



The integral with the induced potential SU{x) from 
Eq. (EH), lu, is 



Cmax kL{6 - {kLf) COS fcL - 3(2 - (kLf) sin kL 
hi = — — 



I2d 



[kLf 



(6.9) 



Now we can perform the summation over modes for 5<t, 
Eq. (|6.6p . analytically in two asymptotic cases: kpL <C 1 
and kpL ^ 1. 

For kpL ^ 1, the evanescent modes give the most 
important contribution to the conductivity^, 



w r 

aikpL « 1) = ^ j 



dx W 

— 72- = r7' (6-10) 

, cosh X t^-Lj 



and to the correction to the conductivity, 

6a{kpL « 1) = ^-I^^ikpLfl , (6.11) 
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FIG. 11: Dependence of the relative correction to conductiv- 
ity on the maximum deformation, Soi^rnax) / for constant 
kpL — 40, in the first order in (,max/d, Eq. (|6.6p . the curve 

2, and the exact transfer matrix solution of integral equa- 
tion, the curve 2. Both the expansion and summation of the 
(7 are not vafid for £,maxkFL/d > 1 for small param- 
eter as mentioned in the text. The dashed line, the curve 

3, is the summation result. The exact solution shows linear 
dependence on deviation even for rather large deviation, and 
this linear dependence is close both to the correction Eq. (|6.6p 
and to the correction averaged over fast oscillations, ^max/'id. 
Parameters of the fiake are^. 



with 

/ = 2 



c?a:sinha;(a;(6 -f- x^) cosha; — 3(2 -I- x'^) sinhx) 
x"* cosh^ X 

0.124. The relative correc- 



and its numerical value is / ; 
tion to conductivity reads 



— = —^C^fL) 

a bd 



0.021- 



'-{kpLY^ 



.12) 



For kpL ^ 1 we average over fast oscillations. In this 
case, only the propagating modes contribute significantly 
to the conductivity. 

To perform the averaging, we replace the summation 
over q by the integration. 



E 



w 

2^ 



dq. 



To simplify subsequent calculations, we make the change 
of variables q = kp sin </>, k = kp cos 0, and then go from 
the integral over dq to the integral over d<j}. The correc- 
tion to the conductivity Eq. (|6.6p has the form 



5(j 



kpW i„ 

ir/2 



COS (/)sin^ (t)siiY{kpL cos < 



(cos-^ (p coa^ {kpL cos ip) + sin (kpLcosip))^ 

In this expression only the term with 3 sin kL / kL 
from Eq. (|6.9p survived: All terms with cos kL vanish 
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after averaging, and the term with —6smkL/{kL)^ is 
smaller than one which is taken into account). For 
large kpL the terms cos{kpL cos(j>) and sm{kpLsm(t)) in 

Eq. (|6.6p oscillate very rapidly. We can represent /J^^^ dcj) 
as a sum of fast oscillating terms, with each term being 
an average over the period, 



with kpLsm(j)n = 2Tm. The integrand / is determined 
by structure of Eq. (|6.6p and Eq. ()6.9p . 



27T f2-K 

f {<!>) = 



sin xdx 



la Jo (a^ cos^ X + sin^ a;)^ 

What is left is the sum over n, 
^max kpW ^ 



2tt 
a 



S(7 



d TT 

Xn = 7r{2n + 1). Finally, 

6a = 



27r^ Vl-(^«/fcFL)2 



'-kpW. 



3.13) 



(6.14) 



The conductivity after averaging over fast oscillations be- 
comes a = kpW/A, and the relative correction to the 
conductivity is 



i+i/2i cos(fci?L COS (j)),shi{kpL cos 4>)), 



5 a ^„ 



2d 



From general physical considerations about the correc- 
tion (see main text), one also expects the dependence 
6a /a ^ £,,nax/d for 6a /a. 



Concerning the correction due to the pseudomagnetic 
fields, it is of the second order in 6 A, and the analytical 
expressions are too cumbersome. Instead, we illustrate 
our conclusions using the numerical solution of the in- 
tegral equation (|6.3p . It is done by multiplying transfer 
matrices for small intervals of the length 6x. Convergence 
with the size of 6x is reached. 
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